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ABSTRACT 

In hierarchical models of structure formation, an early cosmic UV background (UVB) is produced 
by the small (T v ir Si 10 4 K) halos that collapse before reionization. The UVB at energies below 13.6eV 
suppresses the formation of stars or black holes inside small halos, by photo-dissociating their only 
cooling agent, molecular H2. We self-consistently compute the buildup of the early UVB in Press- 
Schechter models, coupled with H2 photo-dissociation both in the intergalactic medium (IGM), and 
inside virialized halos. We find that the intergalactic H 2 has a negligible effect on the UVB, both because 
its initial optical depth is small ( 0.1), and because it is photo-dissociated at an early stage. If the 
UV sources in the first collapsed halos are stars, then their UV flux suppresses further star-formation 
inside small halos. This results in a pause in the buildup of the UVB, and reionization is delayed until 
larger halos (T vlr ij 10 4 K) collapse. If the small halos host mini-quasars with hard spectra extending 
to ~lkeV, then their X-rays balance the effects of the UVB, the negative feedback does not occur, and 
reionization can be caused by the small halos. 

Subject headings: cosmology: theory - early universe - galaxies: formation - molecular processes - radia- 
tive transfer 



1. INTRODUCTION 

Recent progress in high redshift observations has made 
it possible to estimate the global star-formation history of 
the universe to redshifts as high as z ~ 5. Although the 
evolution of the star-formation rate (SFR) is relatively 
well-established at lower redshifts (z <, 2), there is still 
substantial uncertainty about the 2 <, z <, 5 era. A study 
of Lyman break galaxies in the Hubble Deep Field shows 
a decline of the SFR towards high z, indicating a peak in 
the SFR around z ~ 2 (Madau et al. 1996). On the other 
hand, a recent census of z > 3 galaxies, detected in Lya 
emission in a ~ 200 times wider field, suggests that the 
SFR is either constant, or increasing from z ~ 2 to z ~ 5 
(Steidel 1999). At present, these uncertainties render ex- 
trapolations of the observed SFR to still higher (z <; 5) 
redshifts unreliable. 

We do know, however, that some stellar activity must 
have taken place at z > 5. The presence of heavy elements 
observed in the high redshift Lya forest (Tytler, Cowie), as 
well as the reionization of the intergalactic medium (IGM) , 
require sources of metals and ionizing photons in addition 
to the currently known population of galaxies and quasars. 
Such ultra-high redshift objects are indeed expected in 
the current best-fit cosmological cold dark matter (CDM) 
models, which predict that the first virialized halos appear 
at redshifts as high as z ~ 30. It is natural to identify 
these small halos as the sites where the first stars or mas- 
sive black holes are born, forming the first generation of 
" mini-galaxies" or " mini-quasars" . 

The first collapsed halos have masses near the cosmo- 
logical Jeans mass, M ~ 10 4 ~ 5 M Q , implying a virial 
temperature of a few hundred K. In order to form ei- 
ther stars or black holes, the gas in these halos must 



be able to cool efficiently. At this low temperature, the 
only mechanism in the chemically simple primordial gas 
that satisfies this requirement is collisional excitations of 
H2 molecules. At still higher temperatures (T <; 10 4 K) 
and masses (M ~ 10 8 M Q ), efficient cooling is enabled by 
atomic hydrogen. The presence or absence of H2 molecules 
in the first collapsed halos therefore determines whether 
these small halos can form any stars or quasars, before 
larger halos dominate the collapsed baryonic fraction. As 
a consequence, the H2 abundance determines when the 
"dark age" ends: if there is sufficient H2 inside the first 
small halos, the dark age ends earlier; if H2 is absent, it 
ends at the later redshift when the larger scales collapse. 

The fraction of H2 molecules in the post-recombination 
IGM is xh 2 = ^Ha/^H ~ 10~ 6 . In the absence of a back- 
ground radiation field, this fraction rises to ~ 10 -4 at 
the dense central regions of collapsed halos (Haiman et al. 
1996; Tegmark et al. 1997), high enough to satisfy the 
cooling criterion. However, H2 molecules are fragile, and 
are photo-dissociated efficiently by a low-intensity UV ra- 
diation. In a previous paper (Haiman et al. 1997, hereafter 
HRL97), we showed that the flux necessary for H 2 photo- 
dissociation is several orders of magnitude smaller than 
the flux needed to reionize the universe. Based on this 
result, we conjectured that the H2 abundance in the first 
collapsed halos was suppressed soon after a small number 
of these halos formed stars or turned into mini-quasars. 

This negative feedback would imply that the bulk of the 
first mini-galaxies or mini-quasars would reside in dark 
halos with masses at least 10 7-8 Mq, and that the "dark 
age" would end only at the typical collapse redshift of 
these halos. The first condensations at the ~ 1000 times 
smaller Jeans mass would then remain neutral, gravita- 
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tionally confined gas clumps, until they merge into larger 
systems, or are photo-evaporated by the UVB established 
after reionization (Barkana & Loeb 1998). A small fraction 
of these clumps is expected to survive mergers and photo- 
evaporation until z ~ 3; based on their column density 
and abundance, Abel & Mo (1998) have identified this sur- 
vived fraction with the Lyman limit systems observed near 
z ~ 3. A direct observational study of the first collapsed 
halos will be feasible in the future with the Next Gener- 
ation Space Telescope (NGST): with its ~ InJy imaging 
sensitivity, NGST will be able to detect a ~ 10 6 M Q halo 
with an average star formation efficiency at z ~ 10 (Smith 
& Koratkar 1998). 

In this paper, we study the conjecture of HRL97 in a 
cosmological setting, by quantifying the radiative feedback 
of the first collapsed halos on their H2 abundance, as the 
early UVB is established in Press-Schechter models. The 
two main questions we address are: (1) How does the H2 
feedback effect the typical masses and collapse redshifts of 
the first halos; and (2) Could the first collapsed dark halos, 
with masses near the Jeans mass, contribute to the early 
UVB, and to the reionization of the universe ? 

In HRL97, we concluded that the negative feedback sup- 
presses the fragmentation and collapse of small halos. In 
the present paper, we confirm this basic conclusion, al- 
though we find that it depends sensitively on the spectrum 
of the early sources. Our calculations improve HRL97 in 
several ways. We solve the coupled evolution of the forma- 
tion of new halos (using the Press-Schechter formalism), 
the build-up of the early UVB, and the evolution of the 
H2 abundance in the IGM and inside each collapsed halo; 
we explore the sensitivity of our conclusions on the spec- 
tral shape, specifically the X-ray to UV flux ratio; we 
replace our assumption of homogeneous slabs by the cen- 
trally condensed profiles expected for virialized halos; and 
we employ a more realistic star-formation criterion. 

This paper is organized as follows. In § 2, we quantify 
the spectrum of the early UVB after processing through 
the intergalactic H and H2. In § 3, we compute the H2 
abundance in virialized halos under this UVB, and define 
our criterion for star-formation. In § 4, we compute the 
evolution of the UVB coupled with the H2-feedback, and 
present these results in § 5. In § 6, we discuss how these 
results are modified in the presence of early mini-quasars 
with spectra extending to the X-rays. Finally, § 7 summa- 
rizes the main conclusions and implications of this work. 

Throughout this paper, we adopt a background 
ACDM cosmology with a tilted power spectrum, 
(O , Oa, fib, h, o- 8h -i , n) = (0.35, 0.65, 0.04, 0.65, 0.87, 0.96); 
however, our results are relevant in any model where 
the cosmic structure is built hierarchically from the 
"bottom up". For convenience, we also adopt the fol- 
lowing terminology: a "small halo" is defined to be 
a halo with virial temperature T < 10 4 K, or equiva- 
lently, with mass M <, 3 x 10 7 M [(l + z)/ll]~ 3 / 2 . 
Conversely, a "large halo" is defined as a halo with 
virial temperature T 10 4 K, or equivalently, with mass 
M <; 3 x 10 7 M [(l + z)/ll]- 3 / 2 . 



2. RADIATIVE TRANSFER OF THE SOFT UV 
BACKGROUND 



In hierarchical models of structure formation, small ob- 
jects appear at high redshift, and larger objects form later, 
with the total fraction of collapsed baryons monotonically 
increasing with time. The radiation output associated 
with the collapsed halos gradually builds up a cosmic UV 
background. Of interest here is the background in the 
11.18-13.6 eV range, enclosing the Lyman and Werner 
(LW) bands of molecular H2 . Photons in this range do not 
ionize neutral H atoms, and can travel to a large fraction 
of the Hubble length across the IGM, and photo-dissociate 
molecular H2 both in the IGM and inside other collapsed 
halos. Since the presence of H2 is a necessary condition for 
the formation of stars or a central black hole, this leads to 
a feedback that controls the fraction of collapsed baryons 
turned into luminous sources, and the evolution of the 
UVB. Before addressing the feedback itself, in this section 
we quantify the spectral shape of the UVB in the relevant 
frequency range. This is determined by the absorption of 
the flux by intergalactic H and H2 . 

2.1. "Sawtooth" Modulation From Neutral H 

Consider an observer at redshift z b s , measuring the 
UVB at the local frequency v bs, where 11.18eV < ft.f bs < 
13.6eV. The photons collected at the frequency ^ b s are 
arriving from sources located at z s > z Q b s , and have red- 
shifted their original frequency v s > v \ }S to . If during 
its travel, the photon's frequency equals the frequency of 
any atomic Lyman line, it is absorbed by the neutral IGM, 
because the optical depth of the IGM in these lines are ex- 
ceedingly high, e.g. TLy Q ~ 10 6 at z ~ 20. This statement 
applies to all Lyman lines with 2 < n <, 150. The Lya ab- 
sorption (n = 1) has only a small effect on the background, 
since most of the absorbed photons are simply re-emitted 
at the same frequency, except those converted to lower fre- 
quency photons by two-photon decays. In any case, the 
Lya resonance is at 10.2 eV, outside the frequency range 
of interest here. The very high Lyman lines (n <; 150) 
have small enough cross-sections so that the total optical 
depth of the IGM at these frequencies drops below unity. 
However, these resonances have a negligible contribution 
to the total absorption in the 11.8-13.6eV range. Here we 
make the simple assumption that a "dark screen" blocks 
the view of all sources at redshifts above z max . The maxi- 
mum redshift is given by 

1 ^obs ^obs 

where v- x is the frequency of the Lyman line closest from 
above to the observed frequency, v a \, & . 

The processing of the early UVB by the neutral IGM 
leads to a characteristic spectral shape, resembling a "saw- 
tooth" , since the closer v b s is to a resonance, the nearer 
the screen is to the observer. The resulting spectrum is 
shown in a specific model (described below) in Figure [TJ 
The redshift distribution of sources is important for the 
magnitude of this "sawtoothing" effect. If a large fraction 
of the sources contributing to the UVB are located near 
the observer, then the effect is smaller; if the sources are 
spread out over a wider range of redshifts, then the effect 
is larger. In the hierarchical models, new halos form ex- 
ponentially rapidly at high redshifts. At any given time, 
most of the sources are expected to have formed recently, 
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Fig. 1. — The spectrum of the UV background at redshifts z = 30 and z = 15 in the cosmological model of Figure |?| The solid lines show the 
spectrum after absorption by neutral H in the high-redshift IGM; the dashed lines show the spectrum without this absorption. For reference, 
at the bottom of the figure we show the lo catio n of the 76 LW lines of H2. The length of each line indicates its relative contribution to H2 
dissociations (see the Appendix and Figure A16 for the numerical values of the dissociation fractions and oscillator strengths. 



i.e. at redshifts only slightly above that of the observer. 
Since the formation rate of new halos decreases with red- 
shifts, the sawtoothing will be increasingly pronounced at 
later times. This is illustrated in Figure [fl: the flux is sup- 
pressed by a larger average factor at z = 15 than it is at 
z = 30. 



2.2. Modulation From H2 Molecules 

In addition to absorption by neutral H, the UVB may 
also be modulated by the intergalactic H2. The H2 in the 
IGM removes photons from the UVB when these molecules 
are photo-dissociated, and decreases the rate of photo- 
dissociations inside collapsed halos. It is important, there- 
fore, to assess the magnitude of the modulation from the 
intergalactic H 2 . 

When a LW photon is absorbed by an H2 molecule, there 
is on average a ~ 15% chance that it causes a photo- 
dissociation of the molecule. This is the average proba- 
bility that the molecule decays from its excited state to a 
continuum of two distinct H atoms (Solomon, cf. Field et 
al. 1966; Stecher & Williams 1967). We will assume here 
that in the other ~ 85% of absorptions, the LW photons 
is re-emitted at the same frequency, and there is no net 
effect. This is justified because the quadrupole transitions 
among the excited electronic states are much slower than 
decay back to the electronic ground state, and also because 
the small dipole moment that preferred the upper state in 
the absorption will preferentially populate the initial lower 
state in the decay (Herzberg 1950). 

If the H2 fraction in the IGM is xh 2 — "■H 2 / ri H, then the 
total optical depth of the IGM at redshift z and frequency 



v to H2 -dissociations is given by 

2 76 f dt 

T H2 {v, z) = X H2 V" /diss,i/osc,i / dz'c—TlK{z')(t>(Vg>,l^), 

rrieC^ J z dz' 

(2) 

where the sum is over the 76 LW lines of H2; njj = 
8.5 x 10 -6 £lb^ 2 (l + z) 3 is the neutral H number density in 
the IGM; /di 8 s,i is the probability that absorption into the 
i th LW line is followed by a dissociative decay; / osc ,i is the 
oscillator strength of the i th LW line; v z i = v(l+z')/{l+z); 
cdt/dz is the line element in our ACDM cosmology; and 
()>{y, v{) is the absorption line profile centered around the 
resonance frequency v\. We included the equilibrium ortho 
to para ratio of 3:1 of H2 implicitly in equation (Q) by ap- 
propriately multiplying the oscillator strengths with 0.75 
or 0.25 respectively. Furthermore, we assume that the line 
shapes follow a Voigt profile with thermal broadening by 
the IGM temperature (taken to be T = 0.0135(1 + zf K, 
although the total optical depth, and therefore our results, 
are insensitive to this choice) . More information on the de- 
tails of H2 pre-dissociation are given in the Appendix. 

Equation ^ gives a formal sum total optical depth of the 
76 LW resonances, each of which is at a different frequency. 
From the point of view of a photon traveling across the 
IGM, these resonances would be encountered at different 
redshifts. During its lifetime, a UV photon can therefore 
be subjected only to a subset of the LW resonances: as 
we argued in § 2.1 above, the maximum redshift interval 
a UV photon can travel before it is absorbed by a neutral 
H atom corresponds to the redshift between two neutral H 
resonances. A simple calculation of the optical depth us- 
ing equation |^ would result in an overestimate by a factor 
of a few (depending on z and v) of the true value rele- 
vant to the observed modulation of the UVB. To impose 



4 



0.4 



0.3 



m 0.2 



0.1 




Fig. 2. — The total H2 optical depth of the IGM for an observer at redshift z = 30 and z = 15. We assume a constant fixed xn 2 = 2 X 10 -6 , 
and sum the optical depth of the Lyman- Werner lines of H2. At each observed frequency, we only include sources between z ^ s and the nearest 
atomic Lyman resonance. The dashed lines show the effective optical depths averaged in the 11. 18-13. 6eV range. 



the fact that the UVB is simultaneously modulated both 
by the intergalactic H and H2, one needs to change the 
summation in Equation]! from i e(l,76) to i G(i m in,*max), 
where i m ; n is the LW resonance closest in frequency from 
above to is, and i max is the LW resonance farthest above 
v before en count ering an atomic Lyman line (see Figure [I] 
and Figure A16 for visualizations of the relative location 
of the H and - lines). 

In Figure @, we show the resulting optical depth as a 
function of the observed frequency, for two different red- 
shifts, z = 15 and z = 30. As the figure shows, the optical 
depth is largest at ~12.1eV. This is because the largest 
frequency gap between the atomic lines occurs between 
Ly/3 and Ly7, enclosing the largest subset of the H2 LW 
lines with signi fican t contributions to H2 dissociations (see 
Figures g and |A16| ). Fi gure also shows that the maxi- 
mum optical depth is 0.35 and 0.13 for z = 30 and z = 15, 
respectively (assuming ih 2 = 2 x 10 -6 ). Interestingly, 
the characteristic features in Figure @ depend only on the 
relative frequencies of atomic vs. molecular resonant ab- 



sorption lines (shown in Figure A16), i.e. only on atomic 
and molecular physics. 

The total optical depths are near, but somewhat below 
unity. This result has several consequences. It is useful 
here to draw an analogy with the ionization of neutral H. 
The total optical depth of the IGM to H photo-ionizations 
(just above 13.6eV) is r H ~ 100[(1 + z)/21} 3 , which justi- 
fies the simple picture of isolated HII regions (" Stromgren 
spheres"), each bounded by a sharp ionization front, and 
continuously expanding around each ionizing source. The 
situation with H2 is different in two ways: the optical 
depth is less than unity, and is distributed among sev- 
eral LW lines, each of which is individually optically thin, 
and causes dissociations at a different redshift relative to 
each source. As a consequence, no sharp H2 dissociation 
"front" exists. 



A schematic sketch of four distinct regions in the ra- 
dial structure of an H2 photo-dissociation region is shown 
in Figure |^. Region I extends to the usual H I ioniza- 
tion front or Stromgren radius. Inside this region, the 
temperature is high (~ 10 4 K), hydrogen is ionized, H2 is 
photo-dissociated, and the H2 abundance tends to zero. 
Region II is defined to extend from i?; to the radius 
R p ~ 2.5 h' 1 / 2 ^ + z)3o 3/4 A^ 47 kpc, where the H 2 pho- 
todissociation time becomes longer than the Hubble time. 
Here Nlw,47 = A r r J vy/10 47 s _1 is the photon number flux 
in the Lyman Werner bands per second. In this region 
the temperature is that of the IGM, and the H2 fractions 
are still low, but slowly rise to its background value of 
~ 2 x 10~ 6 . Since the optical depth in the Lyman Werner 
lines are very small, the dissociation front is significantly 
diluted. Region III is defined to extend from R p to the ra- 
dius R s « 2/i^ 1 [(l + z)/31]~ 1 - 5 Mpc where Lyman 7 shifts 
to Lyman/3. In this region, the source contributes to the 
soft UV background, but its flux is diminished by absorp- 
tion into the atomic Lyman series. In region IV beyond 
i? s , all photons in the Lyman Werner bands are absorbed 
by an atomic Lyman line. Hence, at such large scales, 
an individual source does not contribute to the UVB and 
intergalactic H 2 dissociation. 

Although the maximum optical depths from the inter- 
galactic H2 are relatively small, the UVB flux can be sup- 
pressed by as much as 1 — exp(— 0.35) ~ 50%. However, 
we are interested in the suppression of the H2 photo- 
dissociation rate when the UVB falls on a collapsed halo. 
This suppression is smaller than 50%, since the photo- 
dissociation rate is given by a weighted average of the ob- 
served flux in the 11. 18-13. 6eV range. We therefore define 
an "effective" optical depth < r > 2 , that directly charac- 
terizes the suppression of the H 2 photo-dissociation rate, 
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Fig. 3. — Schematic sketch of the radial structure of a cosmological photo-dissociation region. Region I is the usual Stromgren sphere, 
bounded by the H I ionization front at Hi. Inside this region, the temperature is high (~ 10 4 K) hydrogen is ionized, and the H2 abundance 
tends to zero. Region II extends from Ri to the radius R p where the H2 photodissociation time becomes longer than the Hubble time. Here the 
temperature is that of the IGM, the H2 fractions are still low, but slowly rise to the background value of ~ 2 X 10 — 6 . Region III extends from 
R p to the radius R B , where Ly7 shifts to Lyman/3. In this region, the source contributes to the soft UV background, but its flux is diminished 
by absorption into the atomic Lyman series. In region IV beyond R B , all photons in the Lyman Werner bands are absorbed by an atomic 
Lyman line. Hence, at such large scales, an individual source does not contribute to the UVB and intergalactic H2 dissociation. See § 2.2 for 
a more detailed discussion. 



as follows: 



fcdiss.uO) = ^XH 2 nH(z)^ 



Ei=il„ /di SS ,i/osc,i J v dv'<j){v', Z^i) to^i 



(3) 



fixed amplitude, and postpone the treatment of the cou- 
pled evolution of the UVB, and the H2 abundance and star 
formation to § 4 below. 



Ei=il, n /di ss ,i/osc,i J v di/(j>(i/, exp[-T H2 (i/, z)], 

(4) 

< T > z = - log[fc diss (z)/fc diss , (2:)]. (5) 

Here J u (z) is the UVB flux before processing by the in- 
tergalactic H2 , and th 2 (v, z) is the optical depth at this 
frequency shown by the solid lines in Figure |^. When the 
sources are assumed to be distributed according to the 
Press-Schechter model (see below), the maximum reduc- 
tion in the H2 photo-dissociation rate is 10% at z = 30 
and 4% at z = 15 (see the dashed lines in Figure ||). 
Note, however, that these numbers assume a constant 
xh 2 = 2 x 10~ 6 . In reality (as will be shown below), the 
intergalactic H2 fraction drops sharply soon after the ap- 
pearance of the first sources, before the UVB builds up to 
the level of causing any H2 feedback inside the collapsed 
clumps. This is because the H2 dissociation time in the 
IGM, 1/fcdiss becomes shorter than the Hubble time at 
J ~ 10 -3 , while the feedback on the dense clumps typi- 
cally turns on only later, when the UVB reaches J ~ 10~ 2 . 

In summary, we have shown that the intergalactic H2 
might cause a maximum of ~ 50% suppression of the 
UVB near 12 eV, but the overall effect on the H2 photo- 
dissociation rate is small: less than 10%, assuming a con- 
stant intergalactic H2 fraction of 2 x 10~ 6 . We will show 
below that the intergalactic H2 fraction is substantially re- 
duced, and therefore this small modulation disappears, by 
the redshifts of interest - when the UVB builds up to the 
levels to photo-dissociate H2 molecules in dense collapsed 
clumps. 



3. H2-COOLING IN THE PRESENCE OF THE UV 
BACKGROUND 

In this section, we consider the behaviour of individ- 
ual, centrally condensed gas clouds, when illuminated by 
an external UVB. Here, we assume the UVB to have a 



3.1. The H2 Abundance in Collapsed Clouds 

The radiative efficiency of a gas cloud by collisional ex- 
citations of H2 molecules depends on the gas temperature, 
density, and H 2 abundance. Although the typical aver- 
age overdensity of a virialized spherical perturbation is 
S ~ 100 — 200, in a centrally condensed cloud, the central 
density can reach a much higher value. This is impor- 
tant in the present context, because the formation of H2 
molecules can be enhanced by the accelerated chemistry 
inside the central, dense regions. In addition, the total in- 
tegrated column density of a centrally condensed cloud is 
higher than that of a homogeneous sphere with the same 
mass and density. Central condensation therefore tends to 
make the clouds more self-shielding against the external 
UVB, and help to preserve the internal H2 molecules. 

To characterize the central condensation, here we adopt 
the density profiles of truncated isothermal spheres (TIS, 
Shapiro et al. 1999). These profiles are especially conve- 
nient, since they have a one-to-one correspondence with 
the final virialized state of the spherical top-hats used in 
the Press-Schechter formalism. The spheres are isother- 
mal at the virial temperature, their average overdensity is 
~ 130, and their central overdensity in the core reaches 
~ 10 4 . The TIS profiles are similar to the spherically av- 
eraged universal profiles found in numerical simulation of 
clusters (NFW, Navarro et al. 1998). The main differ- 
ence from NFW is that the TIS has a central core (the 
existence of cores are indicated by recent simulations of 
galaxy groups, e.g. by Kravtsov et al. 1998). 

We assume that the TIS is composed of the nine species 
H, H~, H + , He, Hc+, Hc++, H 2 , H+, and e~, with the ini- 
tial abundances given by their post -recombination values 
(Anninos and Norman 1996). Our results are insensitive 
to the precise initial conditions, such as the electron frac- 
tion (x e ~ 10 -4 ) and molecule fraction (xh 2 ~ 2 x 10~ 6 ). 
Although these abundances represent the chemical compo- 
sition of the IGM, and we start with a cloud that is already 
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Fig. 4. — The time— evolution of the H2 profile in a truncated isothermal sphere in the presence of an external flux. The example shown is 
for a truncated singular isothermal sphere with T v ; r = 10 3 K, 2 co n = 20, and J21 = 10 — Increasing numbers correspond to equal time-steps 
in logarithmic time, upto a Hubble time of ~ 5 X 10 15 sec. At the outskirts of the cloud the H2 fraction decreases with time. In the core 
molecules are formed initially via H~ + H — > H2 + e~ at a rate exceeding the photo-dissociation rate. At the point where the electron 
fraction drops below 10~ photo-dissociation becomes dominant and H2 is destroyed in the already self— shielding core. 



collapsed and virialized, this should not effect our final re- 
sults. In 1-D spherical simulations (Haiman, Thoul & 
Loeb 1996), we have found that the abundances stay near 
their initial IGM values until the formation of the virial- 
ization shock; soon after virialization the abundances loose 
memory of their initial values. 

We solve the subsequent time evolution of the chem- 
ical abundances, as a function of radius throughout the 
cloud, assuming that the cloud is illuminated by the ex- 
ternal UVB. The UVB flux is assumed to be a power law 
below 13.6eV, 

Jv = -hi ( — ) x sawtooth hv < 13.6eV, (6) 

and to be zero above 13.6eV, since the flux from stellar 
sources would suppressed by several orders of magnitude 
by the neutral IGM. The sawtooth modulation by the in- 
tergalactic H is described above, and the constant ampli- 
tude J21 is in units of 10 -21 erg s _1 cm -2 Hz -1 sr _1 . Our 
results are insensitive to the shape of the UVB, since we 
only utilize the flux between 11.18 and 13.6 eV. The possi- 
bility of a non-stellar flux, extending to higher energies will 
be discussed in section § 6 below. The chemical reaction 
rates and cross-sections are taken from Haiman, Thoul & 
Loeb (1996). We have solved the radiative transfer across 
the cloud, assuming spherical symmetry, and included the 
self-shielding due to both H and H2 . The initial conditions 
are specified by the virial temperature T vlI , and collapse 
redshift z co \\ of the TIS, and the amplitude of the UVB, 
hi. 

An illustrative example of the evolution of the H2 abun- 
dance is shown in Figures ^ and H. Figure [| shows the 
radial profile of xn 2 at 10 different time-steps, equally 
spaced in logarithmic time, during a time-interval that 



corresponds to the Hubble time at z = 20, t w 5 x 10 15 
sec. The cloud is assumed to be a TIS with z co \\ = 20 and 
T V i r = 10 3 K, and the amplitude of the UVB is taken to be 
J21 = 10~ 15 . As Figure [| shows, initially the H2 fraction 
rises in the core: this is because new H2 molecules form 
via the charge transfer reactions H + e~ — > H + 7 and 
H + H — > H2 + e~, utilizing the initial residual free elec- 
trons. By t = 10 14 sec, xh 2 reaches a value near ~ 10 -4 
(cf. the lines labeled 1-6). At this point, however, the 
free electron fraction has dropped by a factor of ~ 2, and 
although the core is self-shielding, photo-dissociation be- 
comes faster than H2 formation. As a consequence ih 2 
drops continuously, to a level below the initial ~ 10~ 6 by 
t = 10 15 sec (cf. the lines labeled 7-10), while x e drops to 
- 2 x 10~ 6 . 

The evolution of xr 2 at the core radius r — tq is also 
shown by the middle solid curves in Figure I for various 
J21. The shell at r = ro is normally representative of 
the whole core, because the density, as well as xh 2 , is flat 
within r < ro (in some cases, however, the H2 fraction may 
still rise within the core, see § 6 below). For comparison, 
in this Figure we also show the evolution of Xb_ 2 m the 
same shell under less intense ( J21 = or 10 -2 ) or stronger 
(J21 = 10~ 4 or 10~ ' 5 ) UV backgrounds. In the absence of 
any UVB, the H2 fraction rises continuously, and reaches 
10 -3 in a Hubble time (top curve); when a flux is turned 
on and increased (bottom 4 solid curves), the H2 fraction 
is continuously suppressed. 

By performing similar calculations for clouds with virial 
temperatures and collapse redshifts in the ranges 10 2 < 
T v ; r < 10 4 K and 5 < 2 co n < 50, we have found that the 
evolution shown in Figures |i| and || remains qualitatively 
the same in all cases. However, the level of flux required 
at which the H2 molecules are photo-dissociation depends 
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Fig. 5. — Evolution of the H2 abundance in the core of the singular isothermal sphere shown in Figure ^; in the presence of various UV 
background fluxes (solid lines). Also shown is the H2 abundance at which the cooling time equals the time elapsed since the collapse of the 
cloud (short dashed line); and the H2 abundance at which the LW lines would become self-shielding (at an integrated H2 column density of 
5 X 10 cm" 2 , long dashed line). For J > 10 — 1,5 the cloud becomes self— shielding, but still cannot cool fast enough. 



on the values of both z co \\ and T v j r . 

3.2. The H2 Abundance and Criterion for Star 
Formation 

In order for a cloud to fragment into stars, or continue 
collapsing to form a central black hole, it is necessary for 
the gas to radiate efficiently. A rough estimate of the 
"critical" H2 fraction can be obtained by requiring that 
the cooling time is shorter than either the dynamical time 
(HRL97), or the Hubble time (Tegmark et al 1997). Here 
we adopt a somewhat more elaborate criterion for star for- 
mation, as follows. For a given combination of z co n, T v i r , 
and J21, the cloud is assumed to be able to cool if, at any 
time during its evolution, the cooling time in the core ex- 
ceeds the "present time" , defined as the time elapsed since 
the formation of the cloud at z C oii- We assume this to be 
both a necessary and a sufficient criterion. 

This definition explicitly takes into account the be- 
haviour of the H2 fraction shown in Figures ^ and ||, es- 
pecially the fact that the H2 fraction in the core first rises 
and then declines again. Our criterion imposes the require- 
ment that the cloud not spend more than a cooling time 
at a given H2 fraction, before the H2 abundance is further 
reduced by photo-dissociation. This definition avoids star 
formation in a situation when the cooling time exceeds the 
dynamical time, but only for a brief interval that would 
be too short for a significant contraction to take place. In 
Figure 0, our criterion is implemented by requiring that 
the solid line (the evolving H2 fraction) goes above the 
short-dashed line (the critical H2 fraction at which t coo \ 
equals the present time t). 

Interestingly, we find that the cloud may undergo a brief 
interval during which it is self shielding in the H2 LW 
bands, but then may become optically thin again as the H2 
abundance drops due to subsequent photo-dissociation. 



This is demonstrated in Figure the long dashed line 
shows the value of xn 2 that would imply an integrated 
H2 column density of Nn 2 = 5 x 10 14 cm~ 2 across the 
cloud, at which the LW lines become self-shielding. As 
the figure shows, if the flux is below J21 ~ 10 -0 7 , the 
H2 fraction does temporarily rise above this critical value; 
but for J21 <; 10 -1 ' 5 , our criterion for star formation is 
still not satisfied. This implies that self-shielding by it- 
self does not guarantee that the molecules are preserved 
and enable efficient cooling. Requiring the cloud simply 
to self-shield would therefore not be a good star forma- 
tion criterion. However, the decline in the H2 fraction is 
due to the decreasing electron abundance which lowers the 
H2 formation rate. Hence, if an additional process keeps 
the e~ fraction greater than ~ 10~ 5 , more H2 could be 
formed. This possibility is discussed in greater detail in 
§6. 

3.3. The Minimum Flux for Negative Feedback 

In addressing the effect of the feedback from the UVB, 
the main question is: given a collapsed cloud with z co \\ 
and T v i r , what is the critical value of the background 
flux, J cr ;t so that star-formation in this collapsed cloud 
is suppressed? In the example shown in Figures |^ and ^ 
(z C oii, T V ir) = (20, 10 3 K), J crit can be read off directly from 
Figure 0, J cr i t = 10 -15 . In order to be able to investigate 
the evolution of the UVB in the Press-Schechtcr models, 
we have repeated the same calculation for a grid of values 
for 5 < z C oii < 50, and 2 < log(T vir /K) < 4. For each pair 
of values, we have found the critical flux for preventing 
star formation, J cr it- 

The resulting J crit = Jcrit (zcoll, T vil ) are shown in Fig- 
ure [| In general, the higher the density (or z co \\) or tem- 
perature, the higher the flux needs to be to prevent star 
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Fig. 6. — The critical flux above which star formation in a halo with virial temperature T v ; r and collapse redshift z co n is prevented. Halos 
smaller than T v ; r as 10 2,4 K cannot cool even in the absence of any flux. Halos larger than T v j r as 10 3 ' 8 K can cool via atomic H, irrespective 
of the H2 abundance. For comparison, the dashed curve show the flux corresponding to one ionizing photon per H atom, the minimum flux 
necessary to reionize the universe. 



formation. We have also found that for temperatures be- 
low log(T V j r /K) < 2.4, the value of the flux becomes irrele- 
vant for star formation, because the cooling function of H2 
drops, and H2 cooling is inefficient at these low tempera- 
tures, even in the absence of any flux. Similarly, the value 
of the flux is irrelevant for log(T v i r /K) > 3.8, because at 
these high temperatures, cooling from (collisional excita- 
tions of) neutral H always dominates over H2 cooling, i.e. 
stars can form irrespective of the H2 abundance. 

The values of J cr it shown in Figure ^ are relatively low. 
For reference, the dashed curve in this figure shows the 
value of J cr it that corresponds to one ionizing photon per 
hydrogen atom, i.e. the minimum flux necessary to reion- 
ize the universe. Since the values of the flux necessary to 
prevent star-formation are 2~4 orders of magnitude below 
this minimum reionizing flux, a negative feedback neces- 
sarily turns on prior to reionization. If the star formation 
efficiency was a universal constant inside all halos that can 
cool efficiently, independent of redshift and halo size, then 
the negative feedback would turn on long before reioniza- 
tion - when the collapsed fraction of baryons f co \\ is ~4 
orders of magnitude smaller then it is at reionization - and 
it would last until / co ii rises by a factor of ~ 100. 

Since new sources appear exponentially fast in the 
Press-Schechter formalism, the relatively "long" interval 
in /coil can translate into a "short" interval in redshift. 
However, it is clear from the above that the fundamental 
"evolutionary parameter" for the negative feedback, fol- 
lowed by reionization, is / co ii, rather than the redshift to 
which it translates to. Our results may therefore vary in 
different cosmologies when stated in terms of redshift, but 
would be relatively robust to change in the background 
cosmology, when stated in terms of / co ii- Finally, we em- 
phasize that the actual redshift (and / co ii) interval be- 
tween the turn-on of the negative feedback and reioniza- 



tion would be shorter if the star-formation efficiency in 
small halos is intrinsically much smaller than in large ha- 
los (see discussion below). 

4. H 2 -FEEDBACK AND THE BUILDUP OF THE COSMIC 
UV BACKGROUND 

In this section, we put together the results of the previ- 
ous sections to compute the redshift evolution of the UV 
background, the intergalactic H2 fraction, and the mass- 
scale of the H2 feedback. The latter is more conveniently 
expressed by defining a critical virial temperature, T CT n,(z), 
that corresponds to the minimum size of a cloud that sat- 
isfies our star formation criterion under a fixed UV back- 
ground J21 and collapse redshift z co n- The evolution of 
T CT i t (z) and the UVB amplitude, J2i(z) are then coupled 
together, since T cr it(z) is directly determined by J2i(z), 
and conversely, J2i(z) depends on the T CT i t (z) above which 
we allow halos to contribute to the UVB. As argued above, 
the evolution of the intergalactic H2 can be considered sep- 
arately, since it does not modulate the UVB at the level 
that would alter the H2 photodissociation rate inside ha- 
los. 

Our model is based on the Press-Schechter formalism. 
We assume that the mass-function of dark halos is given 
by the Press-Schechter formula, N ps (z,M). The mass 
function of collapsed baryonic clouds is obtained from 
N ps (z,M) by the simple replacement M — > fib/^tot-W- 
For the smallest masses, near the Jeans mass (M ~ 
10 4_5 Mq), we apply a correction to take into account the 
non-zero gas pressure in these clouds. This is achieved 
by replacing M — > f^Vl^/Q.^M , where < /b < 1 is the 
fraction of baryons that settle in the potential well of a 
dark halo of mass M by redshift z (Haiman & Loeb 1997). 
We use the mass-virial temperature relation of the TIS 
solutions. Although this relation was derived in a stan- 
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dard CDM cosmology, the scaling of T v ir with redshift is 
changed only negligibly in our ACDM cosmology at the 
high redshifts of interest. 

To parameterize the uncertainties regarding the ionizing 
photon production rate per collapsed baryon, we assume 
that each halo (provided it satisfies the star formation cri- 
terion) has a constant star formation rate for t on years, 
during which it turns a total fraction e of its baryonic mass 
into stars with a Scalo IMF. This parameterization could 
correspond to different physical scenarios. For instance, e 
would remain the same if the star formation rate in each 
halo was doubled, but the escape fraction of ionizing pho- 
tons reduced to 50%. Similarly, a nominal e > 100% could 
correspond to a "top heavy" IMF biased towards massive 
stars relative to Scalo, so that the number of UV photons 
produced per collapsed baryons is increased. We do not 
attempt to summarize all the relevant uncertainties here 
(instead, see Haiman & Loeb 1997), and view t on and e as 
representative parameters in the simplest possible model. 
It is important to note, however, that a combination of 
these parameters is constrained by the reionization red- 
shift 6 z SS 40 (Haiman & Knox 1999), requiring that 
a few ionizing photons per hydrogen atom be produced 
within this redshift interval. 

Within this model, we start from a high redshift (z = 
200), and take small redshift steps towards z=0. At each 
step, the background flux is obtained by summing the flux 
of all sources that are within the imaginary redshift screen 
at z max as discussed in § 2.1, and satisfy the star formation 
criterion of § 3.2, 

hx{z,v) = J dz'c—j v >(z'), (7) 

where j v (z) is the total cmissivity from all luminous halos 
at redshift z and frequency v, 

where F (z , T cr i t [z]) is the total fraction of baryons in 
collapsed halos with virial temperatures above T cr i t . 
Here i on is in units of seconds, and jv(z) is in 
10~ 21 erg s _1 cm~ 3 Hz -1 sr _1 . 

Once this flux is obtained, we find the critical virial tem- 
perature numerically for the new flux J 2 i and redshift, 
using the information contained in Figure ||. 

5. RESULTS 

Our main results are shown for six different combina- 
tions of lifetime and star formation efficiency (t on ,e) in 
Figures case, the figure shows the evolution 

of the ionizing flux 1 at 13.6 eV (Ji on ); the average flux in 
the LW bands, {J T mLw)\ the critical virial temperature 
T cr jt; and the intergalactic H2 fraction. 

Figure shows the results in a model with (t on ,e) = 
(10 7 yr, 10%). The three solid curves in the top panel show 
the evolution of Jlw i n our model with feedback (mid- 
dle curve), assuming that the critical virial temperature is 

1 Note that this flux is not actually visible from a random vantage 
13.6eV across the ionized IGM. 

2 Our model shows that reionization occurs at z < 15. However, we 
a few ionizing photons per H atom (Miralda-Escude ct al. 1999). 



fixed at T vir = 10 2 ' 4 K (upper curve), and at T vir = 10 3 8 K 
(lower curve). The solid line in the bottom panel shows 
the evolution of the actual critical temperature, T cr i t , in 
the coupled solution. These curves demonstrate that ini- 
tially the small halos are allowed to form stars, there is 
no feedback, and the flux builds up as in the constant 
T V ir = 10 24 K case. At redshift z sa 35, the feedback starts 
to be noticeable, as the flux builds up to Jlw ~ 10~ 2 , and 
T cr jt rises above 10 2 4 K. By redshift z — 20, T cr ;t reaches 
10 3 8 K. The dot-dashed lines in the top panel show the 
corresponding evolution of the ionizing flux, J; on ; and the 
dotted line shows the minimum flux required for reioniza- 
tion (corresponding to one photon per hydrogen atom). 
As the figure shows, reionization in this model occurs at 
z ~ 15, and at this redshift the contribution from the small 
halos to Ji on is only ~10% percent 2 . The long dashed 
line in the bottom panel shows the evolution of the inter- 
galactic H2 fraction. As expected, xn 2 drops rapidly to 
a negligible value by z = 35, i.e. before the H 2 -feedback 
turns on. This justifies our neglect of the modulation of 
the UVB by the intergalactic H 2 . The short-dashed line 
shows the virial temperature of the mass scale that is just 
Jeans unstable at each redshift. This shows that the ini- 
tial IGM pressure plays a role only for the smallest halos 
at the highest redshifts. 

The next two figures illustrate the effect of changing 
the total star formation efficiency (or alternatively, the ef- 
ficiency of UV photon production per collapsed baryon). 
Figure H is analogous to Figure 0, but e is lowered to 0.3%, 
the smallest possible value that still produces reionization 
by z — 5. The evolution is qualitatively similar to the 
e = 10% case: the H 2 feedback is active during the red- 
shift interval 13 ^ z 5s 30, before reionization takes place 
at z ~ 5. The star formation in small halos is suppressed, 
and their contribution to the background flux at z = 5 is 
only ~ 4%, less than in the e = 10% case. This, however, is 
mostly an intrinsic effect, rather than a result of the nega- 
tive feedback: as shown by the convergence of the top and 
bottom dashed curves, by the late time when reionization 
occurs, the background ionizing flux is anyway dominated 
by the large halos. The figure also shows that the nega- 
tive feedback has a smaller effect than in the e = 10% case 
(the middle and top dashed lines differ by a smaller fac- 
tor). This is because the background flux builds up more 
slowly, and the feedback is delayed to a later redshift. By 
this later redshift, a larger fraction of the small halos have 
collapsed and formed stars. 

Figure demonstrates the opposite effect of increas- 
ing the star formation efficiency to a nominal e = 100%. 
In this case, reionization occurs earlier, at z 19, and 
the H 2 feedback operates at the earlier epoch between 
24 <; z <> 42. As shown by the large difference between 
the top and middle dashed curves, the negative feedback 
is more pronounced, and it suppresses star formation in a 
larger fraction of the small halos. The total contribution 
of small halos to the ionizing flux at z = 19 in this case 
is ~ 9%. In the absence of the feedback, the small halos 
would dominate the flux by a factor of more than 10 at 
z = 19, and reionization would occur at z ~ 25. 

point until reionization is complete and all sources can see each other at 
assume an equality, because reionization is unlikely to require more than 
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Fig. 7. — Coupled evolution of the background flux J, critical virial temperature T cr ; t , and intergalactic H2 fraction xn 2 . We assume a UV 
production efficiency e=10% and source life time t on = 10 7 yrs. Upper panel: The evolution of the ionizing flux at 13.6eV (dot-dashed lines) 
and the average flux in the LW bands (solid lines). In both cases, the bottom and top curve assume a constant T cr i t = 10 3 and 10 2 4 K, 
respectively, and bracket our actual coupled solution, shown by the thick middle curve. The diagonal dotted line shows the minimum flux for 
reionization, corresponding to one photon per hydrogen atom. Lower panel: The evolution of the critical virial temperature in the coupled 
solution (thick solid curve), and using only the flux from halos with T v ; r > 10 38 K (dotted curve). The long-dashed line shows the evolution of 
the H2 fraction in the IGM. The short— dashed line shows the virial temperature of the mass scale that is just Jeans unstable at each redshift. 




Fig. 8. — Same as Figure^, but assuming e=0.3%. 



Figures |To| — [l2| demonstrate the effect of increasing the 
source lifetime, t on to 10 10 yrs. This long lifetime would 
also mimic the scenario in which constant mergers result 
in multiple recycling of the gas (originally associated with 
a given halo) for star formation. Making the nominal life- 
time longer is equivalent to decreasing the star formation 
rate, and has an effect similar to decreasing e: the buildup 
of the background flux is now delayed (although not even- 



tually suppressed as it is when e is decreased). As a re- 
sult, the reionization redshift is delayed. When e=10%, 
then reionization is at z w 5, and the H2 feedback takes 
place in the interval 10 < z <> 25. The total contribution 
to the reionizing flux from small halos is nevertheless still 
negligible, only about 2%. Figures [ll] and |l2| show what 
happens when e is increased to 100 or 1000% (note that it 
can not be decreased below 10% because then reionization 
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Fig. 9. — Same as Figure but assuming e=100%. 
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Fig. 10. — Same as Figure [?], but assuming t on = 10 10 yrs. 



would occur at z < 5). In both cases, the H2 feedback 
turns on before reionization, and small halos contribute a 
negligible amount to the reionizing flux (5% at z = 10 and 
7% at z = 14, respectively) . 

In summary, Figures [l2] demonstrate that in all cases, 
the H2 feedback suppresses star formation inside the small 
halos, and that the contribution of the small halos to the 
reionizing flux is always below 10%. If reionization oc- 
curs late, then the H 2 feedback makes little difference to 
the overall ionizing flux: the contribution from small ha- 
los to the reionizing flux is then small in any case, since 
by this late time, the collapsed baryon fraction, as well 
as the UV background, is dominated by the large halos 



even in the absence of the H2 feedback. This is the case 
in Figures || and [w| This would be especially true in dif- 
ferent cosmologies, where structure forms earlier, so that 
by the time of a late reionization, the dominance of larger 
halos is even more pronounced. Although not important 
for reionization in these late-reionization models, the H2 
feedback is still expected to occur. On the other hand, in 
models which would predict earlier reionization, primarily 
due to the small halos, the negative H2 feedback is more 
pronounced. In these models, it is the feedback that lim- 
its the contribution of the small halos to reionization, and 
delays reionization until larger halos form. 

Although we have assumed a universal star formation 
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Same as Figure but assuming e=100%. 
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Fig. 12. — Same as Figure |ic|, but assuming e=1000%. 



efficiency, stars in small and large halos are likely to form 
differently, since the cooling mechanisms are different in 
the two cases. Thus, there is no a-priori reason to expect 
the same star formation efficiency in these two different 
environments. There is a preliminary indication from 3-D 
simulations that the efficiency in small halos might indeed 
be small (< 1%, Abel et al). If the star formation effi- 
ciency in the small halos were intrinsically much smaller 
than in the large halos, then the contribution to the flux 
from the small halos would be small, regardless of the Efo- 
feedback. We emphasize that the H2 feedback described 
here would still occur prior to reionization, since the flux 
from the rare large halos alone is sufficient to suppress star 



formation inside the small halos. Note that this scenario 
is different from the one presented in, e.g. Figure |?], since 
in that case the small halos exert a feedback on their own 
population, while here the "feedback" is from a distinct 
population of larger halos. This is demonstrated explicitly 
in Figures 0- 12|, where the dotted curves in each bottom 
panel shows the critical virial temperature T cvli , obtained 
from the UV background due only to the large halos. As 
these figures show, in all cases, T cr it still rises above 10 2 4 K 
and reaches 10 38 K before reionization. 

It is not surprising that the H2 feedback generically oc- 
curs in any model that satisfies the reionization constraint 
by z = 5. This is because, as we have shown in § 3.3, the 
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suppression of the H2 abundance and star formation inside 
clouds with the relevant range of densities and tempera- 
tures requires a UV flux that is only a fraction 10 -2 — 10 -4 
of the minimum flux level required for reionization (one 
photon per hydrogen atom). As the UVB builds up to the 
reionizing level, it must necessarily first reach the level at 
which the H2 feedback sets in. 

6. EARLY MINI-QUASARS AND THE EFFECTS OF X-RAYS 

So far, we have assumed that the early sources have 
"stellar" spectra, and have simply truncated the observed 
flux above 13.6 eV. Since all of our results above depend 
on the flux only in a narrow frequency range below 13.6eV, 
they are indeed insensitive to the precise shape of the spec- 
trum. 

An alternative possibility is that a halo collapsing at 
high redshift forms a central black hole, and turns into 
a " mini-quasar" , with a hard spectrum extending into 
the X-ray regime (Haiman & Loeb 1998). The ionization 
cross-section of neutral hydrogen drops rapidly with fre- 
quency (as i' -3 ), and at frequencies above hv J> IkeV, the 
neutral IGM is optically thin. If the spectra of the early 
sources have a non-negligible contribution at these fre- 
quencies, they would establish an early X-ray background 
(XRB), and change the chemistry qualitatively. 

The importance of the X-rays is that they provide addi- 
tional free electrons, both by direct photo-ionization of H 
and He, as well as indirectly, by collisional ionization of H 
by the fast photo-electrons. The resulting extra free elec- 
trons catalyze the formation of H~ , and therefore tend to 
increase the H2 abundance. In an earlier paper (HRL96), 
we have found that inside dense regions (nn il 1 cm~ 3 ), 
the overall effect of a background flux extending to hard 
X-rays is to enhance, rather than to suppress, the H2 
abundance. The central densities of centrally condensed 
virialized clouds are expected to reach uh ~ 1 cm~ 3 . In 
the case of the TIS solutions adopted here, this density is 
exceeded for halos that collapse at z co \\ 2; 7, and therefore 
a positive feedback is expected when the UV and X-ray 
backgrounds together illuminate these halos. 

In order to illustrate the effect of the X-rays, we have re- 
computed the evolution of the H 2 profiles in the TIS mod- 
els. These computations are identical to the ones discusses 
in § 3.1, except that the flux in equation ^|is modified to 

!J2i (t^) -1 x sawtooth, hv < 13.6eV 

£x j 21 (V)- 1 e (-io 22 [ CTH ,„+o.o8 CTH c„]) i 13 6 _ 10keV 
0, " hv > lOkeV 

(9) 

Here <th,i/ and an c . u are the ionization cross-sections of H 
and He, and the exponential factor mimics the absorption 
above 13.6 eV by the neutral IGM by assuming a fixed hy- 
drogen column density of Ah = 10 22 cm~ 2 and a helium 
column density a factor of 0.08 smaller. The fiducial value 
for Ah was chosen based on a comparison of background 
fluxes computed in our Press-Schcchter model with and 
without absorption by neutral H and He at z = 20. The 
units of J21 are 10 -21 erg s^ 1 cm~ 2 Hz -1 sr _1 . We have 
also introduced a parameter e x to characterize the ratio 
of X-ray to UV flux: for pure stellar sources, e x = 0; for 
typical quasar spectra (Elvis et a!.), e x ~ 1. Values of e x 
in-between and 1 could describe either miniquasars with 



softer spectra, or a mixed population of stars and mini- 
quasars. We have verified that our conclusions are not 
sensitive to the upper energy cutoff, provided this cutoff 
is above ^lkeV. 

In Figure |l3|, analogously to Figure g, we show the evo- 
lution of the H2 fraction at the core radius of a TIS. For 
this illustrative example, we assumed that the TIS has 
z C o\\ — 30 and T vn - = 10 3 K, and that it is illuminated 
by a flux with J21 =0.1 and e x = 0,0.1, and 1, cor- 
responding to increasing amounts of X-ray contribution. 
As in the example shown in Figure ||, in the absence of 
any flux the H2 fraction rises continuously, in this case to 
~ 10~ 3 - 5 (dotted curve). If the UVB is turned on, the 
final H2 fraction is suppressed to ~ 10~ 7 (bottom solid 
curve). However, if increasing amounts of X-rays are also 
added (solid curves, second and third from the bottom), 
then the ^-suppression (relative to the dotted curve) is 
less pronounced. 

We find that a ~ 10% X-ray contribution (e x « 0.1) 
is enough to keep the H2 fraction at a high enough level 
(~ 10 -5 ' 5 ) that satisfies our star formation criterion. This 
is seen in Figure [l3|, as the solid curve labeled e x = 0.1 
reaches the short-dashed curve. For still larger X-ray 
contributions, the final H2 fraction is even higher. It is 
interesting to note that deeper inside the cloud, the sign 
of the feedback is entirely reversed, and the H2 abundance 
is increased relative to the case with no flux. The top 
solid curve in Figure [l^ shows the H2 fraction at the cen- 
ter r = 0, rather than at the core radius, in the e x = 1 
case. Indeed, ih 2 rises above 10 -3 , and reaches a ~ 3 
value than it would in the absence of any background flux. 
Figure |TJ shows the profile of xn 2 across the cloud at vari- 
ous time-steps, and demonstrates that the H2 fraction still 
rises within the core (r < ro) by over an order of magni- 
tude. As described above, the enhancement of the final H2 
fraction is due to the additional free electrons produced by 
the X-rays. The bottom panel in Figure [l^ shows the evo- 
lution of the electron fraction in each case: as expected, 
x e is increased by the presence of X-rays. 

By experimenting with different combinations of z co \\ 
and T v i r , we have found that the situation shown in Fig- 
ure O is typical of collapsed clouds with 5 < z co \\ < 50 and 
10 2 ^K < T vil < 10 4 K. Generally, a - 10% contribution 
by X-rays to the flux is sufficient to keep the H2 fraction 
high enough to allow star formation, and hence to elimi- 
nate the negative feedback caused by the UV background. 
To arrive at this conclusion, we have imposed our star- 
formation criterion at the core radius r = tq. However, we 
have also found that at still smaller radii, the H2 fraction 
is typically much larger, and can rise above the value ex- 
pected in the absence of any background flux. Inside the 
core, net feedback on the H 2 fraction is therefore positive, 
rather than negative. In summary, we expect that if the 
early background had a non-negligible contribution from 
"mini-quasars" with hard spectra, extending to ~ 1 keV, 
then the negative H2 feedback described here would not 
occur. 

7. CONCLUSIONS 

We have examined the build-up of the UV background 
in hierarchical models, based on the Press-Schechter for- 
malism, and its effect on star formation inside small halos 
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Fig. 13. — Effect of X-rays on the evolution of the H2 abundance. The top panel is similar to Figure Q, except that 2 co n = 30, and J21 =0.1, 
and we have added an X— ray contribution to the background flux. The curve labeled J=0 shows the evolution of the H2 abundance in the 
absence of any flux. The label e x on the other curves is a measure of the X— ray/UV flux ratio, where e x = corresponds to a pure stellar 
spectrum with no X— rays, and e x = 1 corresponds to a pure "mini— quasar" spectrum with a spectral index of -1. The bottom panel shows the 
evolution of the free electron fraction in each case. 




Fig. 14. — Effect of X— rays on the profile of the H2 abundance, for the same truncated isothermal sphere as in Figure [13. Increasing numbers 
correspond to equal time-steps in logarithmic time, upto a Hubble time of ~ 3 X 10 15 sec. Note that the H2 fraction stiff rises within the core 
r < rg by an order of magnitude. 



that collapse prior to reionization. Wc have confirmed our 
previous conjecture that there exists a negative feedback 
before reionization. Star formation inside small halos shuts 
off prior to reionization, as an early UV background below 
builds up below 13.6 eV, and suppresses the H2 abundance 
in these halos. 

We have also found that as a result, the contribution to 
the reionizing flux from the small halos is limited to be 



less than a ~10% percent in any hierarchical model. It 
is possible that the star formation efficiency via H2 cool- 
ing in small halos is intrinsically small, as indicated by 
recent 3-D simulations (Abel et al.). The efficiency can 
also be kept low by internal feedback mechanisms inside 
each halo, such as destruction of H2 molecules by an in- 
ternal UV source (Silk 1977, Omukai & Nishi 1999), or 
blow-out of the gas by supernovae (Mac Low & Ferrara 
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1998). In this case, the contribution of the small halos to 
reionization can be negligible regardless of the existence of 
the negative feedback. However, we have shown that the 
negative feedback would still occur from the flux of the 
large halos alone. 

Another possibility is that even though the integrated 
star formation efficiency is high, the star formation rate is 
low, so that the UV flux is built up with a significant delay 
relative to the collapse of the halos. In this scenario, reion- 
ization occurs late (near z ~ 5), when the collapsed baryon 
fraction is already dominated by large halos. Again, in this 
case, the negative feedback would still occur early on, but 
would not have a significant effect on reionization. 

More generally, we have shown that the contribution of 
the small halos to the reionizing flux is small in essentially 
any model where the first collapsed halos do not produce 
hard (E IkeV) X-rays. The reason for this is that the 
flux necessary to suppress the H2 abundance and cooling 
in small halos is 10 2 ~ 4 times lower than the minimum flux 



necessary for reionization, corresponding to an one UV 
photon per baryon. 

Finally, we have shown that the negative feedback would 
be turned into a positive one if there was a significant 
( <; 10%) contribution to the early UV background from 
"mini-quasars" with hard spectra extending to the X-rays 
(upto ~lkeV). In this case, the X-rays catalyze the for- 
mation of H2 molecules by providing additional electrons, 
and the overall effect is to enhance, rather then to sup- 
press, the H2 abundance. The evolution of the small early 
halos is therefore distinctly different in stellar and mini- 
quasar models: this may help shed light on the question of 
what type of sources were responsible for ending the "dark 
age" , and reionized the Universe. 

We thank Jordi Miralda-Escude and Mike Fall for stim- 
ulating discussions. This research was supported by the 
DOE and the NASA grant NAG 5-7092 at Fermilab. Tom 
Abel acknowledges support from NASA grant NAG5-3923. 
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APPENDIX 

THE RELEVANT LYMAN- WERNER LINES OF H 2 MOLECULES 

The bulk of the intergalactic H2 is in the lowest states vibrational and rotational state of the ground electronic state. 
Hence only v=0 with J=0 (para) and 1 (ortho) are populated. There are numerous Lyman Werner band transitions that 
can yield photo-dissociation via a decay from excited electronic states to the continuum of two distinct H atoms. These are 
typically classified as R or P transitions, depending on whether the rotational quantum number changes by -1 or +1 from 
the electronic ground to the electronic excited state (angular momentum must change due to parity conservation) . Hence, 
assuming that both para and ortho hydrogen molecules are abundant in their ground states, all R(0) [i.e. J = — > J = 1], 
R(l) [i.e. J = 1 — > J = 2], and P(l) [i.e. J = 1 — > J = 0] transitions are possible, for each excited vibrational level. 
The above is true both for the Lyman and Werner bands. For the Werner band (X — C 1 n n ), there is an additional 



Q-branch in which A J = is allowedjj. Figure A15 illustrates all possible transitions from the lowest energy states of 
ortho and para H2. 

There is a total of 76 possible transitions from the lowest roto-vibrational state of ortho and para molecular hydrogen 
that are below the Lyman limit of atomic hydrogen (13.6eV). Figure [A16| gives an overview of these 76 lines that are 
relevant for cosmological H2 photo-dissociation. The quantity f osc x fdiss quantifies the contribution of an individual 
line to the overall photo-dissociation. The ./-dependent dissociation fractions in this figure were taken from Abgrall et 
al. (1992). From their work one sees that the dissociation fractions of the Werner Q-branch (which couples to the n - 

3 this is due to the different orbital angular momentum along the inter— nuclear axis of the two states, i.e. E — n for the Werner but £ — E 
for the Lyman transition (A doubling) 
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Fig. A15. — Illustration of all possible absorption lines from the lowest energy states of ortho and para H2. 
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Fig. A16. — The relevant Lyman Werner bands for intergalactic H2 photo-dissociation. The solid columns indicate the value of fdiss x fosc 
and the dotted ones show the dissociation fractions. Lines of the Lyman (Werner) system are labeled by Bi/ (Cv') where v' is the vibrational 
quantum number of the excited state. The location of the atomic hydrogen Lyman series are indicated on the top of the graph. 



state) are extremely small. For the R-branch, however, the fdiss are quite significant (especially for v' = 3 of the Werner 
band). As the figure shows, photons with hv 11.9 eV do not contribute significantly to the dissociation rate due to their 
small f sc- In addition, the Werner bands with v' 3 do not contribute much to the dissociation rate. 



